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ABSTRACT 

In this paper we propose and analyze a distributed algorithm for 
achieving globally optimal decisions, either estimation or detection, 
through a self-synchronization mechanism among linearly coupled 
integrators initialized with local measurements. We model the in- 
teraction among the nodes as a directed graph with weights depen- 
dent on the radio interface and we pose special attention to the effect 
of the propagation delays occurring in the exchange of data among 
sensors, as a function of the network geometry. We derive necessary 
and sufficient conditions for the proposed system to reach a consen- 
sus on globally optimal decision statistics. One of the major results 
proved in this work is that a consensus is achieved for any bounded 
delay condition if and only if the directed graph is quasi-strongly 
connected. We also provide a closed form expression for the global 
consensus, showing that the effect of delays is, in general, to intro- 
duce a bias in the final decision. The closed form expression is also 
useful to modify the consensus mechanism in order to get rid of the 
bias with minimum extra complexity. 

1. INTRODUCTION AND MOTIVATIONS 

Endowing a sensor network with self-organizing capabilities is un- 
doubtedly a useful goal to increase the resilience of the network 
against node failures (or simply switches to sleep mode) and avoid 
potentially dangerous congestion conditions around the sink nodes. 
Decentralizing decisions decreases also the vulnerability of the net- 
work against damages to the sink or control nodes. Distributed com- 
putation over a network has a long history, starting with the pioneer- 
ing work of Tsitsiklis, Bertsekas and Athans 1 1 1 on asynchronous 
agreement problem for discrete-time distributed decision-making sys- 
tems and parallel computing |2|. A simple, yet significant, form of 
in-network distributed computing is achieving a consensus about one 
common observed phenomenon, without the presence of a fusion 
center. Linear average consensus techniques have received great at- 
tention in the recent years |3| — 15|. An excellent tutorial on dis- 
tributed consensus techniques is given in |6|. 

Consensus may be also seen as a form of self-synchronization 
among coupled dynamical systems. In |7][8l, the authors showed 
how to use the self-synchronization capabilities of a set of non- 
linearly coupled first-order dynamical systems to reach the global 
maximum likelihood estimate, assuming reciprocal communication 
links. In particular, in 1 8 1 it was shown that reaching a consensus 
on the state derivative, rather than on the state itself (as in 111 — 1161), 
allows for better resilience against coupling noise. 

The consensus protocols proposed in j3) — (D assume that the 
interactions among the nodes occur instantaneously, i.e., without 
any propagation delay. However, this assumption is not valid for 
large scale networks, where the distances among the nodes are large 
enough to introduce a nonnegligible communication delay. There 
are only a few recent works that studied the consensus problem for 
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time-continuous dynamical systems in the presence of propagation 
delays |i4|, |9| — jT2|. The discrete-time case was addressed in | r|, 
|2, Ch. 7.3], where the authors studied alternative asynchronous lin- 
ear agreement discrete-time algorithms. In particular, |4, 9 | provided 
sufficient conditions for the convergence of a linear consensus proto- 
col in the case of time-invariant homogeneous delays (i.e., equal de- 
lay for all the nodes) and assuming reciprocal communication links. 
The most appealing feature of the dynamical system in 1 4 1 is the con- 
vergence of the state variables to a common value, which is known 
in advance (equal to the weighted average of the initial conditions) 
and delay-independent. Unfortunately, this desired property is paid 
in terms of convergence capabilities, since, in the presence of ho- 
mogeneous delays, the system in |4| is able to reach a consensus if 
and only if the delay is smaller than a given, topology-dependent, 
value. The protocol of |4| was generalized in I10I I12I and 1 11 1 to 
the case of time-invariant nonhomogeneous delays (but reciprocal 
channels) and nonreciprocal channels, respectively. The dynamical 
systems studied in IIOIIIII are guaranteed to reach a consensus for 
any given set of finite propagation delays (provided that the network 
is strongly connected), but this common value is not related to the 
initial conditions of the system by a known function. Similar results, 
under weaker (sufficient) conditions on the (possibly time-varying) 
network topology, were obtained in I2j Ch. 7.3] for the convergence 
of discrete-time asynchronous agreement algorithms. This means 
that the final global consensus achievable by the systems in IIOIIIII 
and |2, Ch. 7.3] is not predictable a priori, so that the protocols in 
the cited works cannot be used to distributively compute prescribed 
functions of the sensors' measurements, like decision tests or global 
parameter estimates. 

Ideally, we would like to have a totally decentralized system 
that reaches a global consensus, for any given set of nonhomoge- 
neous propagation delays (as in 1101 , O Ch. 7.3]) and nonrecipro- 
cal channels (as in [11], |2, Ch. 7.3]), whose final value is a known 
and delay-independent function of the sensors' measurements (as in 
l4l). In this paper we fill this lack and propose a distributed time- 
continuous dynamical system having all the above desired features. 
More specifically, we consider a set of linearly coupled first-order 
dynamical systems, in a network with arbitrary time-invariant topol- 
ogy (not necessarily strongly connected, as opposed to |9| — 112|) 
and nonreciprocal communication channels, modeled as a weighted 
directed graph with weights dependent on the physical radio chan- 
nels. The links among the nodes are affected by time-invariant non- 
homogeneous time offsets, taking into account the propagation de- 
lays, proportional to the relative distances among the nodes, and 
clock offsets among the nodes. Our main contributions are the fol- 
lowing: i) We provide necessary and sufficient conditions ensuring 
local or global convergence of our dynamical system; ii) We derive 
the closed form expression for the consensus, as a function of the 
attenuation coefficient and propagation delay corresponding to each 
link among the sensors; iii) We show how to get a final estimate that 
is not biased by the network geometry and coincides with the glob- 
ally optimal decision statistics that would have been computed by a 
centralized network having a fusion center that has ideal access to all 



the nodes. The most appeahng feature of the proposed system is the 
convergence of the state derivatives to a common known value, for 
any given set of propagation delays and nonreciprocal communica- 
tion channels, with the only requirement that the network be quasi- 
strongly connected. 

2. REACHING CONSENSUS THROUGH 
SELF-SYNCHRONIZATION 

In this section, we first show a class of functions that can be com- 
puted with a distributed approach and then we illustrate the mech- 
anism to carry out the computation without the need of any fusion 
center. 

2.1. Consensus Achievable with a Decentralized Approach 

If we denote by j/i, i — 1, . . . ,N the (scalar) measurement taken 
from node i, in a network composed of A*' nodes, we have shown in 
1 8 1 that it is possible to compute any function of the collected data 
expressible in the form 

^ — 'l— 1 



/(2/i,y2,...,yjv) = h 



(1) 



where {ci} are positive coefficients and {gi} and h are arbitrary 
(possibly nonlinear) real functions on R, and i.e., gi,h : R i-^ R, in 
a totally decentralized way, i.e. without the need of a sink node. In 
the vector observation case, the function may be generalized to the 
vector form 



f(yi,y2 



= h 



(2) 



where {g^} and h are arbitrary (possibly nonlinear) real functions on 
R^, i.e., g^, h : R^ >-* R^, and {C^} are arbitrary square positive 
definite matrices. 

Even though the class of functions expressible in the form (TJ or 
(O is not the most general one, it includes many cases of practical 
interest, like, e.g., the computation of hypothesis testing problem, the 
linear ML estimation fp. TSI, the sufficient statistic in detection of 
Gaussian processes in Gaussian noise, the maximum, the minimum, 
the histograms, the geometric mean of the sensors' measurements, 
and so on |8, 13 1. 

2.2. How to Achieve the Consensus in a Decentralized Way 

The next, most interesting question is how to achieve the aforemen- 
tioned optimal statistics in a totally decentralized network. In this 
paper, we consider a linear interaction model among the nodes, and 
we generalize the approach of 1 8 1 to a network where the propagation 
delays are taken into account and the network is described through a 
weighted directed graph (or digraph, for short), which is a more ap- 
propriate model to capture the nonreciprocity of the communication 
links governing the interaction among the nodes. 

The proposed sensor network is composed of A'^ nodes, each 
equipped with four basic components: i) a transducer that senses 
the physical parameter of interest (e.g., temperature, concentration 
of contaminants, radiation, etc.); ii) a local processing unit that pro- 
cesses the measurement taken by the node; iii) a dynamical system 
whose state evolves according to a first-order differential equation, 
initialized with the local measurements, whose state evolves inter- 
actively with the states of nearby sensors; iv) a radio interface that 
transmits the state of the dynamical system and receives the state 
transmitted by the other nodes, thus ensuring the interaction among 
nearby nodes. 



In the scalar observation case, the dynamical system present 
in node i evolves according to the following functional differential 
equation 

K 

iii't) = giiVi) H ^ "ij — Tij) - Xi{t)) , t>0, 

x,(f) = 0,(t), te[-r,0], i = l,...,N, 

(3) 

where Xi{t) is the state function associated to the i-th sensor; Qiijji) 
is a function of the local observation yi, if is a positive coeffi- 
cient measuring the global coupling strength; Ci is a positive coeffi- 
cient that may be adjusted to achieve the desired consensus; Tij — 
Tij + dij/c is a delay incorporating the propagation delay due to 
traveling the internode distance dij , at the speed of light c, plus a 
possible time offset Tij . The sensors are assumed to be fixed so that 
all the delays are constant. We also assume, realistically, that the 
maximum delay is bounded, with maximum value r — maxij Tij. 
The coefficient Uij measures the amplitude of the signal received 
from node i and transmitted from node j. We assume that the ra- 
dio interface is such that atj = \J Pj\hij\'^ / dl^, where Pj is the 

power of the signal transmitted from node j; hij is a fading coef- 
ficient describing the channel between nodes i and j; is the path 
loss exponent. This requires some form of channel compensation at 
the receiver side, if the coefficients hij 's are complex. Furthermore, 
we assume, realistically, that node i "hears" node j only if the power 
received from j exceeds a given threshold. In such a case, a^j 7^ 0, 
otherwise Uij = 0. The set of nodes that sensor i hears is denoted 
hy Mi — {j = 1, . . . , A'^ : aij 7^ 0}. Observe that, in general, 
Uij 7^ aji, i.e. the channels are non-reciprocal. 

Because of the delays, the state evolution ((Sjl for, let us say, 
t > 0, is uniquely defined provided that the initial state variables 
Xi{t) are specified in the interval from — r to 0. The initial condi- 
tions of (O are assumed to be taken in the set of continuous bounded 
functions (j>i (t) mapping the interval [— r, 0] to R. 

Because of the space limitation, in this paper we focus only on 
the case of scalar observations from the sensors. However, similar 
results can be generalized to the vector case 1131 . 

2.3. Self-Synchronization 

Differently from most papers dealing with average consensus prob- 
lems H] — (6) , where the global consensus was intended to be the 
situation where all dynamical systems reach the same state value, we 
adopt here the alternative definition already introduced in our previ- 
ous work (§]. We define the network synchronization (consensus) 
with respect to the state derivative, rather than to the state. 
Definition 1 Given the dynamical system in a solution {x* {t) } 
o/Q is said to be a synchronized state of the system, if 



k*{t)=u;*, yi = l,2,...,N. (4) 

The system (3j) is said to globally synchronize if there exists a syn- 
chronized state as in <^nd all the state derivatives asymptotically 
converge to this common value, for any given set of initial conditions 
{0J, i.e., 

lim ||x,(f)-tu*|| =0, Vi = l,2,...,Ar, (5) 

tl — VOO 

where 1 1 ■ 1 1 denotes some vector norm and {xi (t) } is a solution of (3}. 
The synchronized state is said to be globally asymptotically stable if 
the system globally synchronizes, in the sense specified in ((j)). The 
system is said to locally synchronize if there exist disjoint subsets 
of the nodes, called clusters, where the nodes in each cluster have 
state derivatives converging, asymptotically, to the same value, for 
any given set of initial conditions {4>i}- 



According to Definition [T] if there exists a globally asymptot- 
ically stable synchronized state, then it must necessarily be unique 
(in the derivative). In the case of local synchronization instead, the 
system may have multiple synchronized clusters, each of them with 
a different synchronized state. In the ensuing sections, we will pro- 
vide necessary and sufficient conditions for the system in {3} to lo- 
cally/globally synchronize, along with the closed form expression of 
the synchronized state. 

3. NECESSARY AND SUFFICIENT CONDITIONS FOR 
SELF-SYNCHRONIZATION 

To derive our main results, we rely on some basic notions of di- 
graphs theory, as briefly recalled next. A digraph is defined as 
<f }, where y is the set of vertices and (f C ^ x ^ is the 
set of edges, with the convention that eij — {vi, vj) G (f if there 
exists an edge from Vj to Vi, i.e., the information flows from Vj to 
Vi. A digraph is weighted if a positive weight, denoted by Uij, is 
associated with each edge dj . The out-degree of a vertex is defined 
as the sum of the weights of all its incoming edges. The in-degree is 
similarly defined. The Laplacian matrix L of a digraph is defined as 
L = D — A, where D is the diagonal matrix of vertex out-degrees 
and A is the adjacency matrix, defined as [A]ij = Uij. A digraph 
is a directed tree if it has A'' vertices and A'^ — 1 edges and there ex- 
ists a root vertex (i.e., a zero out-degree vertex) with directed paths 
to all other vertices. A directed tree is a spanning directed tree of a 
digraph 5f if it has the same vertices of A forest is a collection 
of trees. A digraph is balanced if the out-degree of each vertices is 
equal to its in-degree. A digraph is strongly connected (SC) if any 
ordered pair of distinct nodes can be joined by a directed path. A 
digraph is quasi- strongly connected (QSC) if for every ordered pair 
of nodes Vi and Vj there exists a node r that can reach both Vi and 
Vj by a directed path. A digraph is weakly connected (WC) if any 
ordered pair of distinct nodes can be joined by a path, ignoring the 
orientation of the edges. 

The next theorem is the fundamental result of this paper and it 
provides necessary and sufficient conditions for the proposed decen- 
tralized approach to achieve local/global consensus in the presence 
of propagation delays and nonreciprocal communication links. 

Theorem I (fl3|) Let =\Y , S"} be the digraph associated to the 
network in @, with Laplacian matrix L. Let 'y = [7^, . . . , 7jv]"^ 
be a left ( normalized) eigenvector of L corresponding to the zero 
eigenvalue, i.e., 7"^L = Ojj and \\'y\\ = 1. 

Given the system in (O, assume that the following conditions 
are satisfied: al) The coupling gain K and the coefficients {ci} are 
positive; a2) The propagation delays {rtj} are finite, i.e., Tij < 
T = mSkKi^j Tij < +CXD, Vi 7^ j; a3) The initial conditions are 
taken in the set of continuous bounded functions mapping the inter- 
val [-T, 0] foR^. 

Then, system Q globally synchronizes for any given set of prop- 
agation delays, if and only if the digraph is QSC. The synchronized 
state is given by 



2^. , fiCigiiy^ 



-, Vg, (6) 



li"-i. 



where 7, > if and only if node i can reach all the other nodes of 
the graph through a directed path, otherwise 7^ — 0. 

Theorem [T] has a very broad applicability, as it does not make 
any particular reference to the network topology. If, conversely, the 



topology has a specific structure, then we may have the following 
forms of consensus^ 

Corollary I ( 1131 ) Given system assume that conditions al-a3 
of Theorem\J\are satisfied. Then, 

L The system globally synchronizes to the state derivative 



(7) 



Vg, r = 1, • ■ ■ , A'^, if and only if the digraph contains one 
spanning directed tree, with root node given by node r. 
2. The system globally synchronizes and the synchronized state 
is given by d^J with all 7^ 's positive if and only the digraph 
is SC. The synchronized state becomes 



^ ^ l—l 



V Ci + if y V a. 



(8) 



\/q = 1 , ■ ■ ■ , N, if and only if, in addition, the digraph is 
balanced. 

3. The system locally synchronizes in K disjoint clusters "^i , . . . , 
^ {1, . . . , A'^}£|H'if/! synchronized state derivatives 



E7jCi + A'y^ ^-fiaijTij 



(9) 



Vq G '^k and k = 1, ■ ■ ■ , K, if and only if the digraph '^S is 
WC and it contains a forest with K strongly connected root 
components. 

Remark 1: Robustness with respect to large propagation delays. 

The first important property of the proposed system resulting from 
Theorem [T] is its robustness against propagation delays. It turns 
out in fact that the convergence capability of system ^ is not af- 
fected by the propagation delays. This property represents the major 
difference between our system and the scheme proposed in |4, 9|, 
where, instead, even in the special case of homogeneous delays (i.e., 
Tij = T,'ii ^ j) and undirected connected graph (i.e., aij — aji, 
Vi 7^ j), the average consensus is reached if and only if the common 
delay r is smaller than a topology-dependent threshold, whose value 
is a decreasing function of the maximum graph degree. This implies, 
for example, that networks with hubs (i.e., nodes with very large de- 
grees) that are commonly encountered in scale-free networks, are 
fragile against propagation delays, under the protocol of |4, 9|, even 
in the simple case of homogeneous delays. 

The reason for this difference in the convergence capabilities 
of the two systems in the presence of propagation delays, is a con- 
sequence of the alternative definition of global consensus that we 
proposed for ([3}, with respect to the classical one used in 1 1 1 — 16| 
and |9| — 111|. In fact, we do not require all the state variables to 
converge to a common time-independent value, but only to converge 
towards trajectories given by parallel straight lines. This extra flex- 
ibility provides the additional degrees of freedom that make possi- 
ble the achievement of a consensus on the state derivative without 
requiring any constraint on the propagation delays (besides the ob- 
vious requirement of being bounded). 



'We focus, w.I.o.g., only on WC digraphs. In the case of non WC di- 
graphs, Corollary [T]applies to each disjoint component of the digraph. 

^In general, the clusters "i^i , . . . , ^iffc ^2 not a partition of the set of nodes 
{1,... ,7V}. 



Remark 2: Effect of network topology on consensus structure. 

Theorem [T] generalizes all the previous (only sufficient) conditions 
known in the literature |2 Ch. 7.3], 111, (9)- (l2l for the con- 
vergence of linear agreement protocols in the presence of propaga- 
tion delays. In fact, our theorem provides a full characterization 
of system {3} in terms of necessary and sufficient conditions for ei- 
ther global or local synchronization, valid for any possible degree of 
connectivity in the network (not only for SC digraphs as in |4l, (9) — 
nil ), as detailed next. 

In general, the digraph modeling the interaction among the 
nodes may have one of the following structures: i) f# contains only 
one directed spanning tree, with a single root node, i.e., there exists 
only one node that can reach all the other nodes in the network by a 
directed path; ii) ?f contains more than one directed spanning tree, 
i.e., there exist multiple nodes (possibly all the nodes), strongly con- 
nected to each other, that can reach all the other nodes by a directed 
path; iii) is WC and contains a forest, i.e., there exists no node 
that can reach all the others through a directed path. In the first two 
cases, according to Theorem[T] system ([3} achieves a global consen- 
sus, whereas in the third case the system forms clusters of consensus 
with, in general, different consensus values in each cluster, i.e., it 
synchronizes only locally. In other words, global synchronization is 
possible if and only ;/ there exists at least one node (the root node 
of the spanning directed tree of the digraph) that can send its in- 
formation, directly or indirectly, to all other nodes. In particular, if 
only one node can reach all the others, then the final consensus will 
depend only on the observation taken from that node (see (O). On 
the other extreme, the global consensus contains contributions from 
all the nodes if and only if the graph is SC. If instead no node can 
reach all the others, then the information gathered in each sensor has 
no way to propagate through the whole network and thus a global 
consensus cannot be reached. Still, a local consensus is achievable 
among all the nodes that do influence each other (see 

As an additional remark, the possibility to form clusters of con- 
sensus, rather than a global consensus, depends on the channel coef- 
ficients aij. These may be altered by changing the transmit powers 
Pj. According to the previous comments, the nodes with the highest 
transmit power will be the most influential ones. If, for example, 
we want to write a certain value on each node, we can use the same 
consensus mechanism used in this paper by assigning, for example, 
that value to node i, and use transmit powers such that node i is the 
only node that can reach every other node. 

Remark 3: Closed form expression of the synchronized state. An 

additional important contribution of Theorem[T]is to provide a closed 
form expression of the synchronized state, as given in (|6}, valid for 
any network topology (not only for undirected graphs as in |4|). Ex- 
pression ^ shows a dependence of the synchronized state on the 
network topology and propagation parameters, through the coeffi- 
cients {aij} and the delays {rij}. 

Because of this dependence, the final consensus resulting from ^ 
cannot be made to coincide with the desired decision statistics as 
given by Q, except that in the ideal case where all the delays are 
equal to zero. However, expression (|6f suggests also a method to 
get rid of the bias, as shown in the following algorithm. We let the 
system in ([3} to evolve twice: the first time, the system evolves ac- 
cording to ((Sj and we denote by ijJ*{y) the synchronized state; the 
second time, we set gi{yi) = 1 in ([3} and we let the system evolve 
again, calling the final synchronized state cij(l). From {6](, if we take 
the ratio lo* {y)/uj* (1), we get 



c.*(l) 



(10) 



which coincides with the ideal value achievable in the absence of 
delays. 

Thus, using this simple two-step algorithm, we are able to control 
the value of the synchronized state in advance, without affecting 
the convergence capabilities of the system. This makes our system 
strongly different from the linear agreement protocols proposed in 
Ii2i Ch. 7.3], LlOi .l 1 1. In fact, the dynamical systems studied in ^2) 
Ch. 7.3], IIOIIIII are guaranteed to reach an agreement for any set 
of nonhomogeneous delays (provided that the digraph is QSC in 
Ch. 7.3] and SC in IIOIIIII ). but this common value is not related 
to the sensors' measurements by a known function. In other words, 
the global consensus asymptotically achieved by the protocols in [2] 
Ch. 7.3] and |10, 11 1 is a priori unpredictable, so that the systems 
proposed in the cited works cannot be used directly to distributively 
compute prescribed decision tests or sufficient statistics of sensors' 
measurements, as given in l[T). 

4. NUMERICAL RESULTS AND CONCLUSIONS 

In this section, we illustrate first some examples of consensus, for 
different network topologies. Then, we show an application of the 
proposed technique to an estimation problem, in the presence of ran- 
dom link coefficients. 

Example 1: Different forms of consensus for different topologies 

In Figure [T] we consider two topologies (top row), namely: (a) a 
QSC digraph, (b) a WC (not QSC) digraph with a forest composed 
by two trees. For each digraph, we also sketch its decomposition 
into Strongly Connected Components (SCCs) (each one enclosed in 
a circle), and we denote by RSCC the root SCC of any spanning 
directed tree contained in the digraph. In the bottom row of Figure 
[T] we plot the dynamical evolutions of the state derivatives of sys- 
tem lO versus time, for the two network topologies, together with 
the theoretical asymptotic values predicted by l[6j (dashed line with 
arrows). As proved by Theorem [T] the dynamical system in Figure 
[TJi) achieves a global consensus, since there is a set of nodes (those 
in the RSCC component) able to reach all other nodes. The final 
consensus contains only the contributions of the nodes in the RSCC, 
since no other node belongs to the root SCC of a spanning directed 
tree of the digraph. The system in Figure [IJ5) cannot reach a global 
consensus since there is no node that can reach all the others, but it 
does admit two disjoint clusters, corresponding to the two RSCCs, 
namely RSCCi and RSCC2. The middle lines of Figure [ij)) refer 
to the nodes of the SCC component, not belonging to either RSCCi 
or RSCC2, that are affected by the consensus achieved in the two 
RSCC components, but that cannot affect them. Observe that, in all 
the cases, the state derivatives of the (global or local) clusters con- 
verge to the values predicted by the closed form expression given in 
l|6j-l|9j, depending on the network topology. 

Example 2: Distributed optimal decisions through consensus. 

The behaviors shown in the previous example refer to a given real- 
ization of the topology, with given link coefficients, and of the obser- 
vations. In this example, we report a global parameter representing 
the variance obtained in the estimate of a scalar variable. Each sen- 
sor observes a variable yi — Ai^ + Wi, where Wi is additive zero 
mean Gaussian noise, with variance a^. The goal is to estimate ^. 
The estimate is performed through the interaction system l|3j, with 
functions gi{yi) — yt/At and coefficients Ci = Ai/a^, chosen in 
order to achieve the globally optimal ML estimate. The network is 
composed of 40 nodes, randomly spaced over a square of size D. 
We set the threshold on the amplitude of the minimum useful signal 
so that the underlying digraph be SC. The analog system ^ is im- 
plemented in discrete time, with sampling time Ts = 10~'^ sec. The 





time time 

(a) lb) 
ig. 1. Self-synchronization for two different network topologies: a) QSC 
digraph with three SCCs; b) WC digraph with a two trees forest; Tg = 10""^ 
s,T = 50Ts,K = 30. 

size of the square occupied by the network is chosen in order to have 
a maximum delay t = lOOTs . To simulate a practical scenario, the 
channel coefficients atj are generated as i.i.d. Rayleigh random vari- 
ables, to accommodate for channel fading. Each variable aij has a 
variance depending on the distance dij between nodes i and j, equal 
tcfla?^ =P,/(l + 4). 

In Figure |2l we plot the estimated average state derivative (plus 
and minus the estimation standard deviation), as a function of the 
iteration time. The averages are taken over the nodes, for 100 in- 
dependent realizations of the network, where, in each realization we 
generated a new topology and a new set of channel coefficients and 
noise terms. The results refer to following cases of interest: a) ML 
estimate achieved with a centralized system, with no communication 
errors between nodes and fusion center (dotted lines); b) estimate 
achieved with the proposed method, with no propagation delays, as 
a benchmark term (dashed and dotted lines plus x marks for the av- 
erage value); c) estimate achieved with the proposed method, with 
propagation delays (dashed lines plus stars for the average value); 
d) estimate achieved with the two-step estimation method leading to 

i fTOt (snliri 1 inp*; nlii*; rirrlps fnr the. avpratrp valiip^l Frnm Fi ffure [2l 




time 

Fig. 2. Estimated parameter vs. convergence time, 
we can see that, in the absence of delays, the (decentralized) iterative 



'We use the attenuation factor 1/(1 
undesired event that, for di. 
the transmitted power. 



(fi ) instead of l/(fi- to avoid the 
< 1 the received power might be greater than 



algorithm based on l|3j behaves, asymptotically, as the (centralized) 
globally optimal ML estimator. In the presence of delays, we ob- 
serve a clear bias (dashed lines), due to the large delay values, but 
still with a final estimation variance close to the ML estimator's. In- 
terestingly, if the two-step procedure leading to l llOt provides results 
very close to the optimal ML estimator, with no apparent bias, in 
spite of the large delays and the random channel fading coefficients. 

In conclusion, in this paper we have proposed a totally decentral- 
ized sensor network scheme capable to reach globally optimal deci- 
sion tests through local exchange of information among the nodes, 
in the presence of nonreciprocal communication channels and inho- 
Qogeneous time-invariant propagation delays. Differently from the 
average consensus protocols available in the literature, our system 
globally synchronizes for any set of (finite) propagation delays ;/ 
and only if the underlying digraph is QSC, with a final synchronized 
state that is a known function of the sensor measurements. In gen- 
eral, the synchronized state depends on the propagation parameters, 
such as delays and the communication channels. Nevertheless, ex- 
ploiting our closed form expression for the final consensus values, 
we have shown how to recover an unbiased estimate, for any real- 
ization of delays and channel coefficients, without the need to know 
or estimate these coefficients. If we couple the nice properties men- 
tioned above with the properties reported in [81, where we showed 
that, in the absence of delays, the consensus protocol proposed in 
this paper and in 1 8 1 is also robust against coupling noise, we have, 
overall, a good candidate for a distributed sensor network. 
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